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ABSTRACT 

This paper presents direct evidence for hierarchical galaxy assembly out to redshifts z ~ 3. We iden- 
tify major mergers using the model-independent CAS (concentration, asymmetry, dumpiness) physical 
morphological system on galaxies detected, and photometrically selected, in the WFPC2 and NICMOS 
Hubble Deep Field North. We specifically use the asymmetric distributions of rest-frame optical light 
measured through the asymmetry parameter (A) to determine the fraction of galaxies undergoing major 
mergers as a function of redshift (z), stellar mass (M*), and absolute magnitude (Mb). We find that the 
fraction of galaxies consistent with undergoing a major merger increases with redshift for all galaxies, 
but most significantly, at 5 - 10 a confidence, for the most luminous and massive systems. The highest 
merger fractions we find are 40% - 50% for galaxies with Mb < —21, or M* > 10 10 M at z > 2.5, i.e., 
objects identified as Lyman-break galaxies. Using these results, we model the merger fraction evolution 
in the form: f m (A,M t ,MB, z) = fo x (1 + z) mA . We find tua values ~ 4 — 6 for the most luminous 
and massive galaxies, while lower mass and less luminous galaxies have smaller tua values. We use 
these merger fractions, combined with merger time scales calculated from N-body simulations, to derive 
galaxy merger rates to z ~ 3. We also use stellar masses of HDF-N galaxies to determine the mass 
accretion rate of field galaxies involved in major mergers. We find an average stellar mass accretion rate 
of M G ~ 4 x 10 8 M Gyr -1 galaxy -1 at z ~ 1 for galaxies with stellar masses M* > 10 9 M . This 
accretion rate changes with redshift as: Mq = 1.6 x 10 s (1 + z )° " ±0 - 32 M Gyr -1 galaxy -1 . We also 
find that the fraction of stellar mass density in galaxies involved in major mergers increases with redshift, 
with a peak mass fraction ~ 0.5 for the brightest, Mb < —21, and most massive, M* > 10 10 M , systems 
near z ~ 2.5. By comparing merger fractions predicted in Cold Dark Matter semi-analytic models with 
our results we find a reasonably good agreement for the largest and brightest systems, although we find 
more low-mass galaxy mergers at lower redshifts than what these models predict. 



1. INTRODUCTION 

There are two principle processes in galaxy formation: 
the assembly of mass, both baryonic and dark, through 
accretion and mergers, and the conversion of baryons into 
stars. While the latter process is now statistically mapped 
out to nearly the beginning of the universe using a va- 
riety of techniques (e.g., Lilly et al. 1996; Madau et al. 
1998), we are just beginning to understand the process 
and history of mass assembly. Furthermore, there are sev- 
eral types of mass assembly that are related in an un- 
known way, including the collapse, infall, and accretion 
of dark and baryonic matter, which are likely related to 
how baryonic material is converted into stars. Assembly 
of galaxies and dark halos through mergers and accretion 
is also potentially a major player in black hole and AGN 
evolution, the production of gravitational waves, the trig- 
gering of star formation, and possibly a driver of supernova 
and gamma-ray burst rates. 

Hierarchical mass assembly is also the cornerstone of 
all Cold Dark Matter (CDM) models of galaxy formation 
(e.g., Cole et al. 2000 and references therein). These cur- 
rently favored models clearly predict that dark halos of 
modern galaxies were formed in the past through the pro- 



cess of repeated merging of, and buildup from, smaller 
systems. Although CDM models predict that dark halos 
merge, it is not clear if galaxy formation, or star forma- 
tion, occurs during, before, or after dark-halo mergers. If 
baryons collapse to form stars in dark halos before a signif- 
icant amount of halo merging, then, based on dynamical 
friction arguments, we should witness mergers of galaxy 
stellar components. On the other hand, it is also possi- 
ble that gas cools and forms stars after dark halos merge, 
producing the large galaxies we see today (e.g., Noguchi 
2000). 

One method of determining if and how galaxies form by 
merging is to directly measure the fraction of galaxies un- 
dergoing mergers, and mass assembly occurring by merg- 
ers, at various look-back times and estimate from these 
merger and mass assembly rates. While the star formation 
history of the universe is retrievable in part by examining 
nearby resolved galaxy stellar 'fossil' populations, the mass 
assembly history for nearby galaxies is mostly lost through 
equilibrium and relaxation processes. Some merger tidal 
debris may remain for several Gyrs, such as in our own 
galaxy (e.g., Ibata et al. 2002; Newberg et al. 2002), al- 
though these are likely the result of recent minor mergers. 
There is also considerable observational evidence for re- 
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cent major mergers in the local universe (e.g. Schweizer & 
Scitzer 1988; Borne et al. 2000) and accretion of low-mass 
galaxies onto larger ones (Zaritsky & Rix 1997). However, 
direct evidence for an increase in the galaxy merger rate 
at high redshift has not yet been established, despite the 
considerable circumstantial evidence. This includes evolv- 
ing luminosity functions (e.g., Lilly et al. 1995; Ellis et 
al. 1996), and the appearance of distant 'irregular' galax- 
ies seen in HST WFPC2 images (e.g., Driver et al. 1995; 
Glazebrook et al. 1995; Abraham et al. 1996). Indeed 
an enhanced past merger rate can plausibly explain both 
the " faint- blue-galaxy" excess, and may well be the phys- 
ical mechanism driving the evolution of the blue galaxy 
luminosity function seen at intermediate redshifts. This 
idea was suggested over a decade ago (Broadhurst, Ellis & 
Glazebrook 1992), but direct evidence for evolution in the 
galaxy merger rate has not yet been established. 

The most popular method for measuring the evolution 
of galaxy mergers at high-redshift is through pair counts 
(e.g., Zepf & Koo 1989; Burkey et al. 1994; Patton et al. 
1997; Wu & Keel 1998; Le Fevre et al. 2000) or kinematic 
pairs (Carlberg et al. 2000). After correcting for selection 
effects and biases (e.g. Patton et al. 2002), pair-count 
methods can be used to study major galaxy mergers out 
to z ~ 1. At high redshifts, however, this methods be- 
comes difficult and expensive in telescope time, due to the 
many redshifts needed. There has in fact never been a 
measurement of merger fractions, or merger rates, at red- 
shifts z > 1, although detections of high redshift galaxy 
merging have been claimed (e.g. Neuschaefer et al. 1997). 

Perhaps the best way to understand and characterize the 
merging process is to observe high redshift galaxies and 
determine which are undergoing mergers based on their 
stellar light distributions. We argue in this paper, and in 
Conselice (2003), that this can be done using the observed 
structures of galaxies. At high redshifts, we can do this 
using the Hubble Deep Fields (HDF) (Ferguson, Dickin- 
son & Williams 2000). Conselice et al. (2000a) (hereafter 
CBJ00), Conselice et al. 2000b and Conselice (2003) argue 
that galaxies undergoing major mergers can be identified 
through their large structural asymmetries, one aspect of 
the CAS (Concentration, Asymmetry, dumpiness) mor- 
phological system (Conselice 2003). 

Previously, a color-asymmetry diagram was used by 
Conselice & Bcrshady (1999) to determine the fraction of 
galaxies undergoing mergers in the HDF, based on the 
original WFPC2 images, finding a merger fraction of 40%. 
These previous studies are however potentially biased by 
morphological K-corrections, where a galaxy's appearance 
in the rest-frame UV is not necessarily similar to its rest- 
frame optical morphology. This problem is removed in this 
paper through the use of HDF NICMOS images, where 
the rest-frame optical light of galaxies is sampled out to 
z ~ 2.5. 

By using the CAS system on galaxies found in the Hub- 
ble Deep Field North (Williams et al. 1996), and after 
understanding systematics and biases through simulations, 
we are able to measure the major merger history of galaxies 
out to z <~ 3. We use rest-frame B-band asymmetries with 
other structural and photometric indices such as radii, ab- 
solute magnitude (Mb), and stellar masses (M*) to further 
address the question of how galaxy mergers have evolved 
over time. Using this method we find that out to z ~ 1 the 



fraction of galaxies involved in major mergers increases, 
as has previously been noted (e.g., LeFevre et al. 2000; 
Patton et al. 2002). We argue that the fraction of galax- 
ies undergoing mergers is lower for the fainter and lower 
mass galaxies at higher redshifts. The merger fraction does 
however continue to increase with redshift for the brightest 
and most massive systems. Based on this, we compute and 
quantitatively characterize the merger and mass assembly 
rates of galaxies due to merging, out to z ~ 3. We also 
investigate how the fraction of the total stellar mass den- 
sity in galaxies involved in major mergers changes with 
time, magnitude, and mass, and compare our results to 
predicted values from Cold Dark Matter models. 

This paper is organized as follows: §2 is an explanation 
of the data, including redshift information and how the 
asymmetry parameter is measured for HDF galaxies. In 
§3 we briefly explain the basis for our argument that struc- 
tures of galaxies holds information from which a merger 
origin can be derived, §4 gives the basic results of the 
asymmetry measurements and how they correlate with 
other physical properties, including absolute magnitude, 
color, and stellar mass, to argue that we can determine 
the merger history of galaxies out to z ~ 3. §4 further ex- 
plores the comparison of merger histories with Cold Dark 
Matter (CDM) models and §5 is a summary. The cosmol- 
ogy H = 70 km s" 1 Mpc -1 , fi m = 0.3, and tt x = 0.7 is 
used throughout this paper. 

2. PHOTOMETRIC DATA AND PARAMETERS 
2.1. Imaging Data 

The images we use for our morphological analyses were 
acquired with the Hubble Space Telescope as part of 
the Hubble Deep Field North optical and near infrared 
campaigns. Optical data from the Wide-Field Planetary 
Camera-2 (WFPC2) of the HDF-North (Williams et al. 
1996) are combined with Near-Infrared (NIR) observa- 
tions of the same field taken with the Near Infrared Imag- 
ing Camera and Multi-Object Spectrometer (NICMOS) 
(Dickinson 1998, Dickinson et al. 2000). The optical 
data consist of images in the WFPC2 filters: F300W (U), 
F450W (B), F606W (V) and F814W (I). The NICMOS 
HDF-North images were taken in June 1998 with obser- 
vations in the F110W (J) and F160W (H) near infrared 
bands. We also use a deep K-band image obtained with 
the Kitt Peak Mayall 4m. 

From these images a total of 1212 galaxies in the HDF 
were detected with SExtractor at AB magnitudes brighter 
than J = 27, based on the photometry of Dickinson et al. 
(2000). We compute magnitudes, and radii (§2.3) for these 
galaxies within segmentation maps produced by SExtrac- 
tor (Bcrtin & Arnouts 1996). These galaxies constitute 
the sample used throughout this paper, although we ef- 
fectively use only a fraction of these as we only consider 
galaxies with Mb < —18 to avoid strong selection effects 
and biases (§2.3). 

The shapes, sizes, structures, and morphologies of HDF 
galaxies, as observed in these six pass-bands, are affected 
by instrumental effects, such as point spread functions, in 
addition to the intrinsic morphological K-corrections and 
surface brightness dimming. The angular resolution of the 
NICMOS HDF images is poorer than that of the WFPC2 
data because of the longer wavelength diffraction limit and 
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the pixel undersampling of NICMOS Camera 3. The effec- 
tive PSF FWHM is approximately 0.22" . To minimize sys- 
tematic errors associated with comparing galaxies at dif- 
ferent resolutions (CBJ00), we have convolved the WFPC2 
images to match the NICMOS PSF, and made all morpho- 
logical measurements from these PSF-matched data. Note 
also that the angular diameter distance to galaxies varies 
by less than 40% over the rcdshift range 0.5 < z < 3, and 
by < 20 % for 0.7 < z < 3, where the vast majority of 
the galaxies analyzed in this paper lie. Therefore the res- 
olution of our data in linear physical units (e.g., kpc) is 
nearly constant over the redshift and wavelength ranges 
considered here. 

2.2. Redshifts, Photometry and Stellar Masses 

Photometry is done in all six HDF passbands, and the 
K-band image, for every galaxy, with photometric redshifts 
derived based on observed spectral energy distributions 
(Budavari et al. 2000). We use these photometric red- 
shifts and spectroscopic ones, when available, to compute 
the rest-frame absolute B-band magnitude, Mb, and rest- 
frame Johnson (B— V) color for each galaxy. Rest frame 
magnitudes and colors are computed by fitting the spec- 
tral energy distributions of each galaxy to templates. The 
best fit template is then used to calculate the (B— V) color 
of each galaxy based on interpolation of its observed SED. 
A total of 157 spectroscopic galaxy redshifts were available 
to us from Cohen et al. (2000), Dawson et al. (2001) and 
K. Adelberger/C. Steidel (priv. comm). Figure 1 shows 
the absolute magnitude distribution of all galaxies to our 
magnitude limit of J = 27, with derived (B— V) colors, 
plotted as a function of redshift (z). 

The stellar masses used in this paper come from the 
analysis of Papovich et al. (2001), Papovich (2002) and 
Dickinson et al. (2003). The stellar masses are derived 
using the seven band photometry described above and by 
assuming an initial mass function (IMF), metallicity, and 
a star formation history. The star formation history can 
be modeled, for example, as instantaneous, "bursty", or 
exponentially declining. While there are several different 
possible initial mass functions (IMF) and star-formation 
histories for these galaxies, we use the model results of 
Papovich (2002) where the IMF is Salpeter, the metallic- 
ity solar, and the star formation history of each galaxy is 
monotonic. The stellar masses are also not computed for 
all the galaxies in the HDF, as many are too faint for a 
reliable computation. 

For this particular choice of IMF, metallicity and star 
formation history, the typical uncertainties on the stellar 
masses are about a factor of two (Papovich et al. 2001, 
Papovich 2002). The uncertainties are larger at higher red- 
shifts, where the photometric data does not reach red rest- 
frame wavelengths, and for the bluest galaxies, where vary- 
ing amounts of stellar mass can be hidden by the bright 
light of ongoing star formation. 

2.3. Asymmetries 

2.3.1. Computing Asymmetries 

The main quantitative structural parameter used in this 
paper is the asymmetry (A) index, one property of the 
CAS physical morphology system (Conselice 2003). Asym- 
metry has been used previously as a morphological param- 



eter for nearby (Conselice 1997; CBJ00) and distant galax- 
ies (e.g. Shade et al. 1995; Abraham et al. 1996; Conselice 
& Bershady 1999). The asymmetry computation we use 
is described in detail in CBJ00. The basic computation of 
A involves rotating and subtracting a galaxy image from 
itself, and comparing the summation of the absolute value 
of these residuals to the original galaxy's flux. 

Our computation differs from the prescription set out 
in CBJ00 through our use of the SExtrator dimensionless, 
'Kron radius multiplier' (k r ), times a factor of the semi- 
major axis (a), to define the area within which asymme- 
tries are computed. The SExtractor Kron radius multiplier 
(k r ) is a dimensionless quantity used to scale the ellipse 
whose semi-major (a) and semi-minor (b) axes are defined 
by first moments of the light profile. This is different from 
the prescription of CBJ00 where the Petrosian (1976) n 
radius was used. However, comparisons of the n radius 
used in CBJ00 and the first moment radii for bright HDF 
galaxies reveals that the maximum difference of k r x a 
(hereafter Kron radius) and the Petrosian radius is very 
small, usually a few pixels (~ 0.5") at most, and always 
different by less than 10%. Asymmetries computed using 
the two methods are also indistinguishable for galaxies in 
the HDF. Simulations also show that defining the radius 
in this manner is robust when a galaxy is moved to higher 
redshifts (§2.3.5). 

Asymmetries are computed for all galaxies in each ob- 
served WFPC2 and NICMOS band. The rest-frame B- 
band asymmetries are then derived using the spectroscopic 
or photometric redshift of each galaxy, and then interpo- 
lating rest-frame B-band asymmetries using the observed 
asymmetries. The rest-frame B-band asymmetry is com- 
puted in this way through a linear combination of the 
asymmetry values in the two filters nearest the rest-frame 
B-band wavelength. For galaxies at z > 2.6 where the cen- 
tral rest-frame B-band filter wavelength becomes higher 
than the central wavelength of the H-band, we assume 
A(B) (rost ) = A(H) (0 bs)- 

There are, however, many biases, random, and system- 
atic errors that must be understood in detail, and ac- 
counted for, before we can reliably interpret and use these 
asymmetry values to derive evolution. The general proce- 
dure for measuring asymmetries, as outlined above, relies 
on a robustly defined center, radius, and background area. 
As the Hubble Deep Field North is crowded with galax- 
ies, but not confusion limited, finding suitable background 
areas to measure sky statistics is not trivial. We examine 
this problem, and the centering and radius issues later in 
this paper, although see CBJ00 for a detailed description 
of these problems using nearby galaxies. To address some 
of these issues, we examine a subsample of the 38 brightest 
HDF galaxies, with M B < -18, at redshifts 0.4< z <0.7 
in the HDF to determine how centering routines, differing 
radii and background subtraction methods affect the mea- 
sured asymmetries. We also measure asymmetries within 
radii defined within the SExtractor region, and we inves- 
tigate possible galaxy contamination in §3.3. 

2.3.2. Asymmetry Systematics with Radius 

The problem with choosing a radius to measure asym- 
metries, and other structural parameters, lies in the trade- 
off between galaxy coverage and noise, including contam- 
ination from other galaxies. CBJ00 determined that the 
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larger the radius, the more representative the computed 
asymmetry index is in comparison with a fiducial 'total' 
asymmetry and the better its asymmetry values correlate 
with other physical properties. On the other hand, if radii 
used for computing asymmetries arc large, then more back- 
ground will be present in the computation, which increases 
the noise on the asymmetry measurement. A larger ra- 
dius also increases the chance of contamination from other 
galaxies, and thus requires careful attention to neighboring 
objects. 

While there is no fool-proof way of determining the best 
radius to measure asymmetries, we can constrain this to 
some extent by examining asymmetries of the relatively 
nearby galaxy sample in the HDF computed at various 
factors of the Kron radii. These measurements are plotted 
in Figure 2 as a function of the Kron radii as measured 
by SExtractor. The asymmetry parameter is measured 
for the 38 nearby galaxies in the HDF at five different 
radii: 0.25 k r x a, 0.5 k r x a, 1.0 k r x a, 1.5 k r x a, and 
2.0 k r x a. The average and Ha variations of the asym- 
metries for these 38 galaxies at these radii are: 0.08±0.06, 
0.13±0.11, 0.24±0.26, 0.34±0.31, and 0.44±0.33. 

As for nearby galaxies (e.g., CBJ00), asymmetries in- 
crease at higher radii. Some of the higher asymmetries at 
larger radii are the result of contamination from nearby 
galaxies. We also find that errors on individual asymme- 
try measurements increase when using larger measuring 
radii. To strike a balance between representative radius 
and noise on individual measurements, we use the radii 
0.5 x k r x a to measure asymmetries. This is also the 
radius that matches the apparent radius when estimated 
by eye. Through tests we found this radius to be similar 
to the r(j] = 0.2) radius used to measure asymmetries for 
nearby galaxies, and these radii can be robustly measured 
for galaxies out to redshifts of z~3 (§2.3.5). 

2.3.3. Search Radii 

CBJ00 found that asymmetries do not significantly de- 
pend upon the search size used to find the minimum asym- 
metry. The search size is the pixel size that the asymmetry 
computation code uses to find the minimum asymmetry 
by computing A values at centers that differ by the search 
size (CBJ00) until the minimum asymmetry is found. We 
tested this on the HDF data using several different search 
radii, and always found the same asymmetry, independent 
of search size. For example, when changing the search ra- 
dius from 0.5 to 0.1 pixels, the values of the asymmetries 
only change by 5 A = 0.04 at most. 

2.3.4. Background Measurement Systematics and Effects 
of Correlated Noise 

The method of background removal is a very important 
issue that we address here in some detail. The asymmetry 
computation as outlined in §2.3.1 needs to be corrected for 
background noise, and this is usually done using a patch 
of sky in the image that contains no part of any galaxy. 
This approach is necessary, as an annulus surrounding a 
galaxy may be contaminated by faint outer regions which 
have dimmed below the noise level (e.g., the Tolman (1934) 
effect). 

The background noise level in our images varies from 
place to place in the NICMOS images, and to a lesser 
extent the WFPC2 data, due to variations in quantum 



efficiency and dark current noise over the detector array. 
For this reason, no one "blank spot" used to constrain 
background noise in the asymmetry measurements can be 
considered wholly representative. To determine the de- 
gree to which this may affect the asymmetries, we show 
in Figure 3 the asymmetries of all the 1212 galaxies mea- 
sured using five different background regions throughout 
the HDF images (listed as asymmetry runs 1 through 5). 
These five background regions are at the same physical 
place in each band, and were carefully chosen to span the 
area of the HDF. We only show the resulting asymmetry 
distributions in Figure 3 for the I and H band asymme- 
tries, but these are representative of the B,V and J band 
asymmetries, respectively. 

As can be seen, the I-band asymmetries are relatively 
constant using the different backgrounds, but the H-band 
asymmetries can and do vary. To overcome this problem, 
the effective asymmetry measured in each band is found 
by averaging the asymmetries for the five different back- 
grounds. The la variation in these measurements is added 
in quadrature to the average asymmetry measurement er- 
ror, and this value is used as the asymmetry error. 

We also performed a series of simulations to deter- 
mine the effects of correlated noise on the measurement of 
asymmetries. This was done by creating fake noise maps 
and placing galaxies into them, and then measuring their 
asymmetries. When we correlate the noise by smoothing 
the initial noise map by some filter and then remeasure 
the resulting asymmetries, we find that the variation in 
A is tiny, 5A ~ ±0.03. This is due to the fact that we 
empirically remove the background by using a patch of 
sky, which in principle matches the noise pattern under 
the galaxy being studied. 

2.3.5. Detection and Simulations of Asymmetry 
Variations with Redshift 

A significant issue that must be addressed in any study 
that compares properties of galaxies at different redshifts 
is the fact that measured properties, and the detection 
of galaxies themselves, change solely due to redshift ef- 
fects. The rapidly increasing luminosity distance of galax- 
ies, with the slowly changing angular size distance, pro- 
duces a (1 + z) 4 decline in surface brightness. The result 
is that objects we detect in the HDF at lower redshifts 
might not be detectable at higher redshifts, thereby inval- 
idating some comparisons. Changes in S /N and resolution 
due to redshift can also mask, or mimic, real evolution. We 
addresses these issues using simulations and apply this in- 
formation to correct our asymmetry measurements, and 
to understand our detection completeness. Note that we 
always use interpolated rest-frame B-band values of each 
galaxy, thus we remove all morphological k-corrcctions. 

Simulations were carried out by using the 38 galax- 
ies at 0.4< z <0.7 which have Mb < —18 as described 
in §2.3.1. The images of these objects in their approx- 
imate rest-frame B-band (observed F606W, or l/-band), 
are simulated as they would appear at various redshifts 
from z = 1 to 3 in the redder HDF filters (at z > 3 we can 
no longer sample rest-frame B morphologies). These sim- 
ulations were done by creating the background for each of 
the HDF bands, with the same noise characteristics, then 
randomly placing the simulated galaxies into these back- 
grounds. The galaxies are reduced in resolution, signal 
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to noise, flux and surface brightness, and convolved with 
either the NICMOS or WFPC2 PSF (see Conselice 2003). 

After these galaxies were simulated, we ran the SExtrac- 
tor detection software using the exact same criteria used 
for the original HDF detections (Dickinson et al. 2003b 
in prep). From this, we are able to determine the detec- 
tion completeness at each simulated redshift, and measure 
the asymmetries and radii of the galaxies detected. Do- 
ing this allows us to better understand the systematics 
produced by non-detections and how asymmetry changes 
due to cosmological effects as opposed to real evolution. 
We simulate this local sample of HDF galaxies to redshifts 
z = 1, 1.5, 2, 2.5, 3, 4, 5, and 6. While we do not discuss 
real HDF galaxy asymmetries for objects above z ~ 3 in 
this paper, we include these higher redshift simulations for 
completeness, and to investigate trends with redshift. 

The results for these simulations are shown in Figure 4 
and are quantified in Table 1. First, all of the 38 objects we 
simulate from low— z remain detectable until z ~ 2.5 when 
the completeness begins to drop (Figure 4) . The measured 
values of the asymmetries generally become lower at higher 
redshifts (Figure 4). The average corrections necessary to 
account for these effects are listed in Table 1 and are gen- 
erally low, with differences 5 A = 0.04 - 0.06 for the redshift 
ranges studied in this paper. 

There is also only a very slight difference in retrieved 
asymmetries for faint galaxies with different magnitudes 
in each redshift range. The fainter galaxies at each red- 
shift range are generally affected by noisie more and hence 
the systematics effects (and corrections) are larger, with 
values ~ 0.02 higher, on average, than the brighter sys- 
tems. As this is usually smaller than the random mea- 
surement errors for these faint galaxies, we do not account 
for this small difference. As discussed in CBJ00, where 
similar simulations are done in terms of S/N ratio and res- 
olution, the asymmetry index is not greatly affected by 
the reduced resolution and lower S/N ratios for galaxies 
with Mb < —18. Figure 4 does show however that by us- 
ing a magnitude limit of Mb = —18 we begin to become 
incomplete at redshifts z > 2.5. 

Our radii, defined as 0.5 x k r x a (§2.3.1), also shows a 
slight decrease when the 38 low redshift galaxies are simu- 
lated to higher redshift and then redetected with SExtrac- 
tor (Figure 4b). For the remainder of the paper we use 
these results to apply slights corrections (Table 1) to the 
measured asymmetries and radii for galaxies found in the 
HDF at z > 0.7. 

2.3.6. Final Asymmetry Values and Errors 

The final asymmetry for each galaxy in every band is 
computed by taking the average of the asymmetry values 
computed using the different backgrounds (§2.3.4), and de- 
pending on the redshift, applying a systematic error cor- 
rection as described in §2.3.5. The error of each asymme- 
try measurement is computed by combining the average 
measured error with the RMS of the asymmetry values. 
The average errors on our asymmetry values remain ex- 
tremely low for the brightest galaxies at Mb < — 21 with 
< SA >= 0.04 out to z ~ 3 (Table 2). All of the random 
errors in fact remain rather low, except for the faintest 
galaxies, with Mb > —19, at z > 2.5 where the errors 
approach < 5 A >^ 0.2. 

The RMS variations of computed A values using differ- 



ent sky patches is almost always lower than the computed 
random errors for each asymmetry measurement. The av- 
erage random asymmetry errors are plotted as a function 
of redshift and magnitude in Figure 5 as solid lines, and are 
listed in Table 2. The final corrected asymmetries values 
are plotted as a function of redshift in Figure 5. 

3. METHODOLOGY 
3.1. The Merger Criterion 

A major goal in contemporary astrophysics is determin- 
ing how galaxies formed and evolved. As we approach this 
goal, it is fair to ask if the structures of galaxies give any 
clues towards solving this problem. Is it possible that the 
morphological appearance of a galaxy is only representa- 
tive of temporary 'weather', and does not relate to the 
fundamental underlying evolution? While we only briefly 
address these questions here, see Conselice (2003) for a 
detailed discussion and an introduction to the CAS (con- 
centration, asymmetry, dumpiness) classification system 
where it is quantitatively argued that galaxy structures 
reveal fundamental information. 

In general, a galaxy's appearance or morphology is de- 
termined by a variety of different effects. Some of these are 
produced by projection, most notably in the case of edge- 
on spiral galaxies, but as argued in Conselice (2003) mor- 
phology is largely the result of physical processes, such as 
star formation, interactions and mergers with other galax- 
ies, and the past history of these events. 

As described in CBJ00 and Bcrshady, Jangren & 
Conselice (2000), by using a computationally consistent 
method to compute asymmetries there are strong and 
physically meaningful relationships with other parameters, 
such as color and concentration (CBJ00). This requires 
using a method for computing asymmetries that is not 
affected by choice of center, radius, or to first order, reso- 
lution. 

To determine how useful asymmetry and other param- 
eters are for distinguishing galaxies in various phases of 
evolution, CBJ00, Conselice et al. (2000b) and Conselice 
(2003) computed asymmetries for over 200 nearby galaxies 
in all phases of evolution, including high-z analogs such as 
starbursts and ULIRGs. When asymmetries are combined 
with color, or other structural information, galaxy types 
such as ellipticals, spirals, and irregular galaxies can be 
roughly distinguished from each other. As expected, el- 
lipticals are symmetric, red objects; while later-types are 
both bluer and more asymmetric (CBJ00). We also ar- 
gue in these papers that galaxies undergoing major merg- 
ers can be distinguished from those evolving quiescently 
through their global asymmetries. 

There are also fairly strong correlations between (B— V) 
colors and dumpiness (S) values with the asymmetry in- 
dex for non-mergers, such that more asymmetric galaxies 
are bluer and have higher dumpiness values. For these 
non-mergers there exists a small distribution of A values 
at at all S and (B— V) values, with a natural scatter a(A)'m 
asymmetries (see Figures 7 & 8 in Conselice 2003). We de- 
fine these distributions in Conselice (2003) and use them 
to identify statistical outliers with high asymmetries which 
we identify as mergers. If we set a limit of A merger > A(S, 
B— V) + 3a (A) for major mergers, then we find A(B) mcrgcr 
~ 0.35 for both the bluest and most clumpy galaxies yet 
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observed. Redder or less clumpy galaxies have a 3a(A) 
deviation less than A = 0.35, thus we are being conser- 
vative with this limit. From a nearby sample of <~ 240 
galaxies, nearly all objects that deviate more than 3a(A) 
from the asymmetry-color and asymmetry-clumpiness re- 
lationships (Consclicc 2003) are galaxies involved in major 
mergers. As such, we use A(B) mcrgor = 0.35 as our limit 
for identifying major mergers. 

Other evidence that ^4(B) mcrgor > 0.35 includes galaxies 
with asymmetries larger than A(B) mcrger showing kine- 
matic evidence for merging based on broadened HI line 
profiles (Conselice et al. 2000b). The measured asym- 
metries of simulated galaxies in major merger simulations 
also have A > 0.35 when undergoing merging (Conselice 
& Mihos 2003 in preparation). These simulated galaxies 
have lower asymmetries with A < J 4(B) merger before, and 
after, merging events. 

3.2. Eddington Bias 

We argue in §4 that galaxies at higher rcdshifts are more 
asymmetric in their rest-frame B-band morphologies. We 
use this to further argue that, at the very least, the most 
massive galaxies must be forming through the merger pro- 
cess. An important question to ask however is how this 
result is biased by our observational random errors, which 
in noisy data can mimic higher merger fractions, analogous 
to the aberration in star or galaxy counts due to random 
errors (Eddington 1913). 

The basic idea behind the Eddington bias is that within 
an intrinsic distribution of some observed quantity, in the 
presence of more and more noise, there will be larger mea- 
surement tails. In star counts the effect is to scatter more 
counts into bright bins from fainter bins. In this section 
we investigate if the increase in A values can be explained 
by increased noise scattering intrinsically low asymmetry 
values into the high asymmetry bins. 

To investigate the importance of Eddington bias in pro- 
ducing higher measured asymmetries for galaxies at higher 
rcdshifts, we plot in Figure 6 the asymmetry values for our 
sample divided into different redshift and absolute magni- 
tudes. We also list in Table 2 the average random error of 
the asymmetry values for the galaxies plotted in each of 
the rcdshift/absolute-magnitude bins shown in Figure 6. 
From these values, and the other information in Figure 6, 
we argue that the Eddington bias is not a major effect 
when comparing asymmetries for galaxies at similar abso- 
lute magnitudes at different redshifts. The errors remain 
below 6A = 0.10 up until z = 3 when they become large. 

To show this, and to understand the limitations of our 
data and the effects of random errors, we carry out Monte 
Carlo simulations, the results of which are plotted on Fig- 
ure 6, as the number in the lower right of each panel. This 
number is the sigma likelihood that an increase in claimed 
asymmetries between two bins is not due to an increase in 
random errors at higher redshifts. That is, it is the sig- 
nificance that any increase in asymmetries is real and not 
due to Eddington bias. 

These simulations are done by assuming that the ran- 
dom errors are Gaussian distributed, with a full width at 
half maximum given by the increase in the error listed in 
Table 2. To determine how significant our observed in- 
crease in asymmetries are at a given magnitude, we add 
in resulting random errors to the asymmetries for galax- 



ies in the lowest redshift bin (0 < z < 1). We then re- 
compute the mean and standard deviation of the resulting 
asymmetry distribution. This allows us to determine the 
likelihood that an increase in average asymmetries, and la 
distributions, at different redshifts are due to an increase 
in random errors induced by being at higher redshifts. 

The result of these simulations is that in almost all cases, 
the statistical probability of a claimed asymmetry increase 
being due to a random configuration of increased errors is 
very small. Note that we do not claim an increase between 
each redshift interval for every magnitude (see §4.4). From 
these simulations we are confident at the 3 - 10 a level that 
all claims for increases in asymmetries between different 
rcdshifts is a real effect, and not a result of Eddington 
bias. The only exception is the change in asymmetries 
seen for the — 20 > Mb > — 21 galaxies between < z < 1 
and 1 < z < 2. This is the only interval in which a claimed 
increase in asymmetries has a significant less than 3 a. 

Note that we do not use the z > 3 bin in our analysis 
and only show it here for comparison purposes. We lat- 
ter quantify the Eddington bias in another way by always 
plotting the ± random error merger fractions. Doing this, 
we find that the trend of increasing asymmetries at higher 
rcdshifts is still present, verifying the statistical argument 
made in this section. 

3.3. Contamination from Projection 

One potential problem with the asymmetry methodol- 
ogy is that occasionally two galaxies will overlap to pro- 
duce a two dimensional image that looks asymmetric. The 
two galaxies themselves can look very smooth and sym- 
metric, but when projected on or near each other they 
could potentially be identified as a major merger through 
a resulting high asymmetry. Although we only measure 
asymmetries for galaxies within well-defined radii pro- 
duced through SExtractor, which deblends objects, there 
are still possible cases of near neighbors that could produce 
a higher signal. 

As such, we must be able to constrain the importance 
of this effect. Overlapping galaxies in the Hubble Deep 
Field are, however, rare, with only a handful of obvious 
cases (White, Keel & Conselice 2000). Through a visual 
examination of all objects identified as a merger through 
asymmetries, only seven are galaxies that might be in pairs 
(Figure 7). We say 'might' here as it is often the case 
that these pairs are real physical associations, and some 
of the photometric redshifts suggest that they are in fact 
two nearby galaxies, perhaps in the early phases of a ma- 
jor merger. In any case, since our sample contains ~ 70 
galaxies identified as mergers, only 10% of our mergers 
have something that could resemble a pair, although many 
of these show evidence for being real physical units due to 
similar colors, and in at least one case similar spectroscopic 
redshifts, and thus are not chance projections. 

4. RESULTS 

4.1. Merger Candidates in the Hubble Deep Field 

All galaxies brighter than Mb = — 20 between < z < 3 
found in the HDF are shown in Figure 7 in the F814W 
band, divided by redshift intervals: < z < 1, 1 < z < 2, 
and 2 < z < 3. The mergers brighter than Mb = — 20 
are also listed in Table 3. The images shown in Figure 7 
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are the appearance of each galaxy in F814W, with three 
numbers overplotted, which are from top to bottom: abso- 
lute magnitude, rest-frame B-band asymmetry (A(B) res t), 
and redshift. Galaxies which are statistically likely to be 
merger candidates, with A(B) > 0.35, have a solid box in 
the upper right corner of their panel. 

This figure shows, among other things, that the asym- 
metry index is able to pick out systems that would be 
chosen as mergers through visual estimates. Systems with 
high asymmetries are however clearly in various phases 
of merging. Some objects, particularly at high rcdshifts, 
appear to be two galaxies beginning to undergo a merger. 
Some also appear to have a low signal to noise morphology. 
This demonstrates the power of the asymmetry index to 
remove ambiguity and subjectivity in determining which 
galaxies are merging. While any given handful of galaxy 
morphologists would pick out different mergers by eye from 
this list, we are free from this concern as our method is 
purely automated, and can be understood statistically. 

As discussed in Conselice (2003), CAS parameters at 
high redshift are best used in an ensemble sense, such as 
finding the fraction undergoing major mergers, while any 
one measurement can be dominated by random errors. For 
example, the asymmetry values listed on Figure 7 some- 
times have random errors as large as 8 A = 0.2, as do galax- 
ies with A(B) < 0.35. These errors can remove, or add, 
galaxies into the merger bin, a fact that we account for in 
§3.2 and §4.4 when analyzing the merger history of field 
galaxies. 

4.2. Physical Properties of Asymmetry and Size 

4.2.1. Asymmetries 

In Figure 5a we plot the rest-frame asymmetries of our 
sample galaxies, brighter than Mb < —18, as a function of 
redshift and luminosity. In Figure 5b we plot the asymme- 
tries as a function of redshift and stellar mass for systems 
with M* > 10 8 Mq . The symbols on Figure 5 represent 
the magnitude or stellar mass of each galaxy, with larger 
symbols representing brighter or more massive galaxies. 
Note that, particularly at high redshift (z > 1), some of 
the galaxies consistent with mergers are relatively bright 
and massive. 

The average asymmetries of objects brighter than the 
given magnitude limit listed in Table 4 are plotted in Fig- 
ure 8 as a function of redshift, where the respective bright 
magnitude limits are labeled next to each line. The av- 
erage asymmetries generally increase from the z = — 1 
range to the z = 1 — 2 range, and then decreases at higher 
rcdshifts. This is especially true at the fainter magnitude 
limits. At brighter magnitudes the average asymmetry 
generally increases with redshift, peaking at z ~ 2.5, and 
declining thereafter. At the brightest magnitude limit of 
Mb < — 22 the average asymmetry increases greatly be- 
tween z ~ 1 and 1.5 and then stays large at higher red- 
shifts. 

4.2.2. Colors 

We can get some idea of the stellar populations that 
make up galaxies seen in the HDF by examining their 
rest-frame Johnson (B— V) colors. These colors are de- 
rived from an interpolation of the broad-band photometry, 
based on the known spectroscopic redshift, or computed 



photometric redshift, for each object. A plot of rest-frame 
color versus redshift is shown in Figure 9 out to z ~ 2.5, 
the limit where we can measure rest-frame (B— V) colors. 
An obvious feature of this plot is the presence of a large 
population of faint and very blue galaxies at redshift above 
z ~ 1.5 which are not seen at lower redshifts. Some of 
these galaxies are systems with low asymmetries seen in 
the corresponding diagram of asymmetry versus redshift. 

Many of these very blue and faint galaxies are however 
the result of cither noisy photometry, imprecise photo-zs, 
bad K-corrections, or a combination of these effects, as 
very few nearby galaxies have negative (B— V) colors. Pop- 
ulation synthesis models also show that objects with mean 
weighted single stellar population burst ages of 50 Myrs, 
with no subsequent star formation, have a color (B— V) 
~ 0.04 (Bertelli et al. 1994). In fact, since most of these 
colors are derived from spectral energy distribution fits to 
broad brand photometry, there are likely to be some sys- 
tematic errors in individual colors. We therefore only use 
these colors in later figures as representative ensembles in 
different populations. The values of these colors also do 
not affect any of the quantitative results of this paper. 
There does appear, based on Figure 9, to be an absence of 
bright red galaxies with Mb < —20 and (B— V) > 0.5 at 
z > 1.5, although we cannot out rule that some intrinsicly 
red galaxies appear blue on Figure 9 due to systematic 
effects. 

4.2.3. Asymmetry and Color Evolution in the Bright 
Galaxy Population 

The relationship between asymmetry, color and mag- 
nitude for HDF galaxies is shown by plotting rest-frame 
B-band asymmetries as a function of Mb into two different 
redshift bins in Figure 10. Figure 10a and Figure 10b show 
Mb vs. A(B) where the points are colored according to 
their rest-frame (B— V) values, within the redshift ranges 
< z < 1.5 (Figure 10a) and 1.5 < z < 3 (Figure 10b), 
respectively. This roughly divides the sample into galax- 
ies dominated by visible star-formation and those that are 
not. Figure 10a shows that the brightest galaxies in the 
low-redshift bin have low asymmetries and red (B— V) col- 
ors, consistent with smooth galaxies with old stellar popu- 
lations, such as ellipticals. There is also a clear bifurcation 
in Figure 10a such that asymmetric galaxies are generally 
blue, and the low-asymmetry objects are generally red (see 
§4.3), the bluer objects possibly the result of merger in- 
duced starbursts. We see the opposite however at higher 
redshifts, 1.5 < z < 3, (Figure 10b) where the bright- 
est galaxies tend to be blue objects undergoing starbursts, 
which have high asymmetries, suggesting they are poten- 
tially major mergers. 

4.2.4. Size Evolution 

Another property that we investigate is the size evolu- 
tion of galaxies in the HDF. For this we use the 0.5 X \Zy X CI 
radius within which asymmetries are measured. Like all 
the other parameters discussed in this paper, the appar- 
ent sizes of objects depends on redshift, as dictated by 
the angular size distance. However, the sizes of objects 
may also appear smaller due to surface brightness dimming 
that makes the outer parts of galaxies difficult to detect. 
To address this issue we investigate, through the simula- 
tions discussed in §2.3.5, how measured sizes of galaxies 
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change due to cosmological surface brightness dimming. 
The change in galaxy sizes due to this effect arc plotted in 
Figure 4b, where average radius differences and their la- 
variations are plotted as a function of rcdshift. 

When we apply this correction to the measured sizes in 
the rest frame B-band, and examine the resulting distribu- 
tion as a function of redshift, we get Figure 11. There are 
two interesting properties in this figure. The first is that 
there appears to be a lower limit on galaxy sizes, which is 
around 3 kpc. This limit is likely the result of the SExtrac- 
tor detection method and effects from the PSF, as galaxies 
smaller than this certainly exist in the local universe (e.g., 
Conselice et al. 2002). There is in fact a bias in the way 
that SExtractor measures the sizes of galaxies, which typ- 
ically depends upon luminosity. Therefore, the best way 
to view Figure 11 is as a relative change in galaxy sizes 
with redshift. The range and distribution of galaxy sizes 
appears roughly constant between 1.4jzj3, while for zjl.4, 
there are are significantly more galaxies with radii above 
7 kpc. Short of a significant amount of cosmic variance, 
this evolution is real, as these sizes have been corrected for 
cosmological dimming. There therefore appears to be no 
galaxies at z > 1.4 larger than 10 kpc in the HDF-N, even 
after correcting for redshift effects. 

This is consistent with the idea, but does not prove, 
that nearby large galaxies are forming from the mergers of 
lower-mass, and presumably smaller, galaxies. An alterna- 
tive interpretation is that these galaxies are still forming 
(inside out) from accreted intergalactic gas cooling into 
stars after z <~ 1.5. Size and luminosity evolution (Ellis et 
al. 1996) alone do not prove that mergers are occurring 
in the galaxy population, but both effects are consistent 
with this idea. For the remainder of the paper we examine 
the structures of galaxies in the Hubble Deep Field and 
argue that mergers are increasingly common out to z <~ 3, 
especially for the most massive systems. 

4.3. Star- Formation versus Merger Formation 

We can do a general test to determine if, and approx- 
imately how much, star formation is the cause of asym- 
metries in HDF galaxies by examining their star forming 
properties as traced by color. If star-formation is respon- 
sible for producing asymmetries of our sample, then we 
would expect more asymmetric galaxies to be dominated 
by star formation. 

We can use the color-asymmetry diagram to investigate 
this question. The color-asymmetry diagram, as intro- 
duced in Conselice (1997), and discussed in CBJ00 and 
Conselice et al. (2000b) is a diagnostic tool that plots 
the disturbance of a galaxy with a measure of its spec- 
tral shape, which signifies the ages of its stellar popula- 
tions. Figure 12 shows the HDF rest-frame (B— V) color- 
rest frame B-band asymmetry A(B rcst ) diagram plotted in 
bins of absolute magnitude, Mb- The differences and sim- 
ilarities between the local galaxy population (sec CBJ00) 
and the HDF galaxies can be seen. The diagonal line is 
the relationship between asymmetry and color character- 
ized by CBJ00 for nearby normal galaxies. The solid part 
is an extrapolation of this relationship to colors bluer than 
(B— V) <~ 0.4, which is typically among the bluest col- 
ors found for e.g., dwarf irregulars in the nearby universe. 
While there are many galaxies with modest asymmetries 
and blue colors, there are few analogs of nearby ellipticals 



in the HDF with red colors and low asymmetries (CBJ00). 

Figure 12 can be used to argue that high asymmetries 
are not likely solely produced by massive amounts of star 
formation. At each magnitude and redshift range there ex- 
ist a range of A rest values for any given color. In fact, the 
galaxies consistent with mergers (those to the right of the 
vertical dashed line) do not have colors that significantly 
differ from blue non-mergers. They are, on average, bluer 
than the total population, but there are always galaxies 
with similar, or bluer, colors that do not have similarly 
high asymmetries (this can also been seen in Figure 10). 

If vast amounts of star formation are responsible for 
large asymmetries, then the galaxies with the bluest colors 
should have the largest asymmetries, which is clearly not 
the case. The galaxies with high asymmetries are unique 
in terms of their structure, which we interpret for the rea- 
sons in §3 and Conselice (2003) that they are undergoing 
major mergers. 

4.3.1. Galaxy Populations in the Color- A symmetry 

Diagram 

As mentioned earlier, there appears to be a lack of cor- 
responding modern day ellipticals, or very red-symmetric 
objects at any redshift bin in the HDF. This may simply 
be due to the fact that ellipticals are bluer at high redshift 
and/or due to the very small co- moving volume of the HDF 
at low redshift which is 320 h~ 3 Mpc 3 at 0.1 < z < 0.5, 
where few, if any, galaxies with L>L* are present (Dick- 
inson et al. 2003). Morphologically selected elliptical- like 
galaxies however do exist in the HDF, particularly around 
z = 1 (Stanford et al. 2003). They are however bluer than 
nearby ellipticals as their stellar populations are younger 
then those in nearby ellipticals (Menantcau et al. 1999), 
an effect also seen in high redshift clusters (e.g. Dickinson 
1997; Stanford, Eisenhardt & Dickinson 1998). Passive 
evolution from z ~ 1 to will create a reddening of about 
0.1 magnitudes, which would place these symmetric bluer 
galaxies in the general area of nearby ellipticals (CBJ00). 

Figure 12 further shows that many galaxies in the 
HDF are potentially undergoing mergers with asymme- 
tries ^4(B)> 0.35. These objects have asymmetries incon- 
sistent with being normal, late-type galaxies undergoing 
quiescent star formation as compared to nearby galaxies. 
As we go to higher redshifts the fraction of galaxies with 
asymmetries consistent with merging increases (Figure 12; 
§4.4). 

What are the galaxies in the HDF that are not identified 
as mergers? At low redshifts, these are simply the mor- 
phologically familiar disks and ellipticals (Figure 7). At 
higher redshifts, the galaxies inconsistent with merging are 
not readily identifiable, through eye-ball estimates, with 
any specific local morphological type, but span a range, 
from very compact galaxies with possible 'tidal' features 
to diffuse blob like systems (Figure 7) (see also Giavalisco, 
Steidel & Macchetto 1996). There are also symmetric sys- 
tems at redshifts z > 1.5, although these are generally 
blue. 

4.4. Galaxy Merger Fractions and their Evolution 

In this section we use the asymmetries of HDF galaxies 
to measure the evolution of implied major merger frac- 
tions out to z ~ 3. As in the previous sections we avoid 



9 



morphological K-corrections by using the rest-frame B- 
band asymmetries of galaxies in the HDF and applying 
the corrections as needed at higher redshifts as described 
in §2.3.5. 

As previously described, we define a major merger as a 
galaxy whose rest-frame B-band asymmetry is larger than 
^merger = 0.35 for the reasons discussed in §3. By taking 
the ratio of galaxies with asymmetries A > A morgor to the 
total number of galaxies in a given parameter range, the 
implied merger fractions out to z ~ 3 can be computed as 
a function of absolute magnitude Mb, stellar mass (M*), 
and redshift (z). This allows us to determine how dif- 
ferent galaxy types have evolved as a function of mass 
and time. These merger fractions are listed and plotted 
in four different lower galaxy absolute magnitude limits, 
M B = -18, -19, -20, -21 in Table 5 and Figure 13, and 
mass limits M* = 10 8 M , 10 9 M , 10 95 M Q , and 10 10 
M Q in Table 6 and Figure 14. Each merger fraction is com- 
puted for galaxies brighter, or more massive, than these 
limits. For example, the Mb = — 18 bin contains all the 
galaxies in the Mb = —19, —20 and —21 bins. 

The solid large circles in Figures 13 and 14 are the in- 
ferred merger fractions based on the asymmetry measure- 
ments using the various magnitude and mass limits shown 
on the panels in Figures 13 and 14. The high and low val- 
ued green crosses are the computed merger fractions after 
adding and subtracting, respectively, the la error from 
each asymmetry measurement, and then recalculating the 
merger fractions. This is a graphical representation of a 
likely outcome of the effect of random errors and system- 
atic biases on the measured distribution of asymmetries 
(i.e., the Eddington bias) (§3.2). 

Merger fractions from the magnitude selected pair stud- 
ies of Patton et al. (1997) and Le Fevre et al. (2000) 
and the kinematic pairs study of Carlberg et al. (2000) 
are also plotted on Figure 13 at their most representative 
magnitude regime. These authors compute merger frac- 
tions by finding the number of galaxies at each redshift 
range separated by some projected distance (usually < 20 
kpc), and in the case of Carlberg et al. (2000), a relative 
velocity difference (< 500 km s _1 ). These authors however 
use various magnitude limits, and it is straight forward to 
compare each of these points to our values. Carlberg et al. 
(2000) use a magnitude limit of Mb ~ —20.5 while Patton 
et al. (1997) and LeFevre et al. (2000) find pairs within 
magnitude limits of Mb = —18 and Mb = —19, respec- 
tively. None of these studies use limiting stellar masses as 
we do here. 

Note that we may be missing a substantial number of 
dusty merging galaxies, such as sub-mm sources, which 
we are not considering because they are too faint to be de- 
tected within our limits. These would increase the merger 
fractions above those calculated here for galaxies within 
our stellar mass limits (Burgarella et al. 2003). They 
would not however effect the derived merger fractions for 
systems within the given magnitude limits. 

4.4.1. Fitting Functions 

We characterize the evolution of these merger fractions 
by fitting a simple power-law increase with redshift: f m (A, 
M*, M B , z) = fb x (1 + z) mA , where f is the merger frac- 
tion at z — and tua is the slope of the merger fraction 
evolution, such that more steeply rising merger fractions 



have higher ttia values. We perform these fits to quanti- 
tatively parameterize how merger fractions evolve, and to 
compare with previous work and theoretical studies that 
also use this parameterization. We have fit this model to 
the data in two different ways: a simple unweighted least- 
squares fit to all the merger fractions and by fitting the 
merger fractions by holding the z ~ point to the value 
found by Patton et al. (1997). 

The results of these various fits are plotted as two dif- 
ferent colored lines in Figures 13 and 14, with the fitted 
parameters listed in Tables 7 & 8. A fit to only the asym- 
metry merger fractions are shown as the straight dark blue 
lines on Figures 13 and 14. The tua values for these fits 
are fairly low for the fainter Mb > — 20 (Table 7) and low 
stellar mass systems with M* < 10 9 5 M (Table 8) with 
typical mA values ~ 0.5 - 1. However we find that for 
the brightest and highest mass systems with Mb < —21 
and M* > 10 10 M Q the merger slopes are quite steep with 
values of 3.7±0.3 and 5.9±1.3, respectively. As we argued 
in §2 and §3.2 this increase in asymmetry between dif- 
ferent redshifts cannot be accounted for by morphological 
k-corrections, detection incompleteness, systematic errors, 
random errors or Eddington bias. 

If we use the Patton et al. (1997) merger results as a 
fiducial z ~ bench-march for all magnitude and mass 
cuts, and then fit the merger fraction evolution based on 
this, we get higher fitted slopes with values near rriA ~ 
1.5 — 2. These fits are shown in Figures 13 and 14 as the 
cyan colored lines. 

We can also examine merger fraction evolution when 
only fitting up to certain redshifts. Tables 7 and 8 list 
the values of m,A for these fits out to z ~ 1 and z ~ 2 
as a function of limiting stellar mass and magnitude. At 
these lower redshift limits the merger fraction slope, tua, 
becomes quite steep, except for the brightest and most 
massive systems (Tables 7 & 8). The slopes of these fits 
are typically tua = 2.5 — 3 for galaxies with Mb > —20 or 
M* < 10 9 ' 5 M© . At z ~ 2 the merger fractions are lower 
and the fitted slopes, tua, decrease to tua ~ 2 for systems 
with M B > -20 or M* < 10 9 5 M Q . We see the opposite 
effect however for the most massive systems, with M* > 
10 10 M© . For these galaxies the fitted power-law slope, 
m,A, is quite low, between z ~ to 1, with a value of tua 
= 1.7. The brighter systems with Mb < —21 also have 
lower mA slopes between z ~ — 1 with values tua ~ 1-5. 

This bifurcation in mA values (Tables 7 & 8) be- 
tween bright/massive and faint/lower-mass galaxies sug- 
gests that massive galaxies form from mergers much ear- 
lier than lower mass systems. The steep decline in merger 
fractions suggests the most massive galaxies underwent a 
massive merger phase in the distant past and have be- 
come the quiescent massive galaxies (or quiescent massive 
galaxy components) we see in the nearby universe. 

4.5. Merger and Stellar Mass Assembly Rates 

Using further assumptions we can investigate the merger 
and mass accretion rates of galaxies from z = out to 
z ~ 3, or back to when the universe was only ^2.1 Gyrs 
old. It should be kept in mind however that the calcula- 
tions we perform below are somewhat speculative as we 
do not yet have a firm understanding of many of the prop- 
erties we are assuming throughout the following analysis. 
One possible problem with this is that the Hubble Deep 
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Field occupies only a small volume of space and a large 
cosmic variance may make these results inapplicable glob- 
ally. We know that in dense areas such as clusters, the 
merging properties are probably different than in the field 
(e.g., van Dokkum et al. 1999). Despite this, the Hubble 
Deep Field samples the high-redshift universe and thus 
within the limitations imposed by the uncertainties in our 
assumptions, we can compute certain evolutionary quan- 
tities for the first time. 

4.5.1. Merger Rates 

To understand the merger and mass accretion rates we 
must know the dynamical time-scale, and mass ratios, of 
mergers which will result in systems with high asymme- 
tries. That is, we have to have some idea of what type 
of major merger will produce a galaxy structure with 
an asymmetry larger than A morgor . The two most im- 
portant parameters for understanding this are: the time 
scale which a galaxy will remain asymmetric, such that 
A > ^4 m crgcr, and the mass ratios necessary to produce 
A > A morgcr . The time-scale for a merger also varies with 
rcdshift solely due to the different physical conditions of 
galaxies in the past, although we ignore this for the mo- 
ment and assume that mergers occur within the same time- 
scale at all redshifts. 

Estimates of mass merger ratios and time scales can be 
computed through N-body simulations of galaxies under- 
going the merging process. These simulations have been 
done and are fully reported in Conselice & Mihos (2003 in 
preparation; hereafter CM03). We briefly summarize the 
conclusions of this study and use this information to de- 
termine the merger and mass assembly rates for galaxies 
seen in the HDF. 

The merger of two galaxies of nearly equal mass remains 
asymmetric for roughly 0.9 Gyrs, that is A > 0.35 during 
this length of time (CM03). There is however a variation 
in this time range, partially for mergers between systems 
that do not contain similar masses. Also, we have no a 
priori method of determining the mass ratios of the sys- 
tems that produced the galaxies which we see as mergers. 
To simplify this, we assume that each merger is produced 
from two galaxies of approximately equal stellar masses. 

Using this modeled length of time as the average time 
scale for which A > ^4 m0 rger, we can compute the merger 
rate as a function of rcdshift. We measure this by assum- 
ing, as earlier, that each galaxy with Aq > ^4 mcr ger is cur- 
rently undergoing a major merger that lasts for ~ 0.9 Gyr, 
within the co-moving volume between each redshift inter- 
val. Using these number, we can then calculate the merger 
rate </>m , defined as the number of mergers occurring per 
co- moving volume (in Gpc 3 ), divided by the time-scale of 
the merger (0.9 Gyrs). 

These merger rates, 4>m, for galaxies brighter than 
Mb = — 19 are plotted as a function of rcdshift (z) in Fig- 
ure 15 as a solid line. We calculate the merger rate as 4>m 
= 3.6 x 10 5 Mergers Gpc" 3 Gyr" 1 between z ~ 0.4 - 0.8 
and find a peak rate of 4>m = 2 x 10 6 Mergers Gpc -3 
Gyr -1 between z = 0.8 and 1.4. At redshifts higher than 
z ~ 1 the merger rate declines towards the merger rate 
values found at z < 1. These rates, as well as the mass 
accretion rates (§4.5.2), are listed in Table 9. From Fig- 
ure 15 the merger rate appears to be nearly constant, but 



the merger fraction generally increases with rcdshift (Fig- 
ure 13-14). This difference is due to the fact that we are 
sampling more volume at higher redshifts, and thus the 
number of mergers per unit volume remains relatively flat. 

4.5.2. Mass Accretion Rates 

The stellar mass accretion rate per co-moving volume 
(M) is plotted as the dashed line and on the right axis of 
Figure 15. The stellar mass accretion rate is computed 
in a similar way to the merger rate, 4>m- This is found 
by assuming that each major merger with A > A merger 
consisted previously of two galaxies of equal mass that 
merged. Therefore the mass accretion rate, M, is the sum 
of half the total mass involved in mergers in a redshift in- 
terval divided by the co-moving volume in that interval, 
divided by the major merger time scale. The resulting M, 
roughly follows the form of the merger rate with a peak 
value of 6.0xl0 6 M Q Gpc -3 year -1 between z — 0.8 and 
1.4. 

It may seem surprising that the mass accretion rate and 
merger rate track each other very well out to z ~ 2.5 
since the masses of high redshift galaxies are on average 
smaller than those at low rcdshift (Dickinson et al. 2003). 
If galaxies involved in mergers were drawn from a random 
distribution, then for a merger rate which is constant with 
z would correspond to a mass accretion rate that declines 
with z. However, because the galaxies involved in mergers 
are of higher mass at higher redshifts, the mass accretion 
rate does not decline. 

The stellar mass accretion rate per galaxy M G is also 
plotted as a function of redshift in Figure 16. The val- 
ues of M G are calculated by determining the amount of 
stellar mass added to galaxies from major mergers in a 
rcdshift interval divided by the total number of galaxies 
within that interval per merger time scale (0.9 Gyr). It 
is the average stellar mass accreted onto a galaxy due to 
major mergers per unit time. 

The stellar mass accretion rate per galaxy, Mq, for 
galaxies more massive than 10 9 , 10 95 and 10 10 M are 
plotted on Figure 16. For galaxies at M* > 10 9 the peak 
rate is 5.5x10 s M Galaxy" 1 Gyr" 1 at z ~ 2.5. The 
peak rate is 1.16 xlO 9 M Q Galaxy -1 Gyr -1 and 4.72 xlO 9 
M Q Galaxy -1 Gyr" 1 for systems at M* > 10 9 5 M and 
M* > 10 10 M Q , respectively. The evolution of M G can be 
further parameterized as: 

M G (z) = 1.6xlO 8 M Galaxy -1 Gyr -1 (l+z)°- 99±0 - 32 (1) 

for systems with M* > 10 9 to z = 2.5. The best fit for 
systems at M* > 10 9 5 is 



M G (z) = 1.8 x 10 8 M Galaxy -1 Gyr -1 (l + z) 1 - 47±0 ' 25 and 

(2) 

M G (z) = 5.3xlO 6 M Galaxy -1 Gyr -1 (l + z) 5 - 3±al6 (3) 

for galaxies with M* > 10 10 . These best fits are plotted 
as dashed lines on Figure 16. 
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4.5.3. Building Galaxies Through Mergers 

Using the derived values in §4.5.2 we can investigate 
how much of the stellar mass of modern galaxies is formed 
through major mergers. This can be done by integrating 
equations 1-3 from z ~ 3 to z ~ where we observe galax- 
ies in their most evolved form. By integrating equations 
1-3 we find that the addition of mass for these galaxies in 
the Hubble Deep Field North due to mergers is at most 
~ 10 10 M Q since z <~ 3. This shows that on average the 
most massive galaxies only slightly double in stellar mass 
due to major mergers. This is not enough stellar mass to 
produce an L* galaxy at z ~ 0, which has a stellar mass 
of ~ 10 11 Mq . However, the vast amounts of ongoing 
star formation in these systems, probably induced by this 
merging, will form stars along with this addition of already 
existing stellar mass obtained through major mergers. The 
HDF may also be devoid of the most massive, > L* , galax- 
ies, as more massive systems at z > 2 do exist (Shapley et 
al. 2001). This is certainly a volume effect and these mas- 
sive systems are rare, and the more typical galaxies found 
in the HDF are those that dominate the mass density. 

Equations (1-3) reveal the first mapping of galaxy stellar 
mass assembly from mergers as a function of stellar mass. 
Combining this with star formation histories, and the re- 
sulting stellar mass formation histories (Dickinson et al. 
2003), it is possible to determine the formation histories 
of galaxies from z <~ 3 as a function of their initial stellar 
mass. 

4.5.4. Mass and Luminosity in Mergers 

We use the above information to determine the fraction 
of stellar mass, and rest-frame B-band luminosity, involved 
in major mergers, within our stellar mass and magnitude 
limits, as a function of redshift. These fractions are calcu- 
lated by determining the total amount of mass and light in 
all galaxies within a given redshift interval, and comparing 
this to the amount of light and mass within that interval 
that are attached to galaxies involved in major mergers. 

Figure 17 shows that roughly 5% to 25% of the stellar 
mass and luminosity in HDF galaxies are involved in ma- 
jor mergers for systems with Mb > —20 or M* < 10 9 5 
Mq . For the highest mass and brightest systems with 
M* > 10 10 M or M B > -21 there is a high mass frac- 
tion peak of 0.5, and luminosity fraction peak of ~ 0.6, at 
z ~ 2.5 and a rapid decline at lower rcdshifts. 

This again demonstrates that high mass and luminous 
galaxies underwent major mergers at high redshift, but are 
not doing so in the nearby universe. Figure 17 also shows 
that although the fraction of stellar mass in the brightest 
and most massive systems involved in major mergers de- 
clines at lower redshift, the fraction of the luminosity com- 
ing from major mergers is relatively constant until z <~ 1, 
when it drops. It should be kept in mind however, that 
the stellar masses used to create Figure 17 are from the 
stars that dominate the light of these galaxies and there 
is a possibility that these galaxies are hiding a significant 
amount of mass in the form of old stars (e.g., Papovich 
2002; Dickinson et al. 2003). These results are however 
consistent with these massive and bright galaxies being el- 
lipticals, or bulges formed through mergers ~ 9 Gyrs ago. 
It is also consistent with peculiar galaxies at higher red- 
shift transforming into normal Hubble types (e.g., Driver 



et al. 1998; Brinchmann & Ellis 2000). 

4.5.5. Possible Biases 

The evolution in Figures 16 and 17 appears steepest for 
objects above the largest mass threshold, i.e. > 10 10 M . 
Given the stellar mass evolution inferred for the general 
galaxy population from Papovich et al. (2001) and Dick- 
inson et al. (2003), it is worth noting that a fixed mass 
threshold such as > 10 10 Mq corresponds to increasingly 
rare and unusual objects at higher redshifts. For example, 
at z <~ 0, 10 10 Mq is ~ 0.1 of the mass of an L* galaxy, 
while at z ~ 3 it appears to be more or less the M* mass 
of Lyman-break galaxies, that is more massive objects are 
evidently rare. Thus, by using a fixed mass threshold, we 
are examining rates for objects that are increasingly ex- 
treme at higher redshift, and which (in hierarchical clus- 
tering theory) are themselves ever more strongly biased (in 
terms of their clustering) relative to the underlying dark 
matter density distribution. 

This might have some implications for the results we 
see in Figures 16 and 17. The steepness of the apparent 
evolution at high masses might be an artifact due to this 
effect, although this is not likely the case. Unlike the mass 
function, the rest-frame optical luminosity function does 
not evolve dramatically with redshift, and yet there still is 
a rapid evolution for objects with Mb < —21 (Figure 16). 
These bright objects, unlike galaxies with M* > 10 10 M Q , 
are not increasingly rare at higher redshift. 

4.6. Comparison to Models 

Many papers have recently made merger history pre- 
dictions based on the CDM paradigm (Khochfar & Burk- 
ert 2001; Gottlober, Klypin & Kravtsov 2001) which we 
compare to our results. The idea that galaxies build up 
from mergers began with the modeling of this process by 
Press & Schechter (1974) and White & Ress (1978). These 
early models traced the abundances of dark halos as a 
function of time with some assumptions for how galaxies 
form through gas cooling and feedback. Later, these mod- 
els were expanded to include an initial power-spectrum 
of dark halos assuming that dark matter is cold (White & 
Frenk 1991). This Cold Dark Matter (CDM) approach us- 
ing an extended Press-Schechter formalism (Bower 1991; 
Bond et al. 1991) has proven to be successful at repro- 
ducing many nearby and high-redshift galaxy properties, 
through semi-analytic and N-body modeling (e.g., Kauff- 
mann, White & Guiderdoni 1993; Cole et al. 1994; Guider- 
doni et al. 1998; Kauffmann et al. 1999; Sommerville & 
Primack 1999; Cole et al. 2000; Somerville, Primack & 
Faber 2001). However, the basic idea that galaxies and 
dark halos merge, the fundamental idea of CDM models, 
has until now never been tested observationally at high 
redshifts where most galaxy formation occurs. 

There are several results from this paper which we can 
quantitatively compare with CDM predictions of structure 
formation. A basic test is to examine how merger fractions 
change as a function of redshift in CDM models, within 
magnitude and stellar mass limits. Comparisons with the 
semi-analytic model results of Benson et al. (2002) are 
shown in Figures 13 and 14 as solid red lines. These pre- 
dicted merger fractions were computed by finding all ma- 
jor mergers, defined, for purposes of this model, as systems 
merging with a mass ratio of 1:3 or higher and within the 



12 



past Gyr, at a given redshift, in the semi-analytic simu- 
lation of Benson et al. (2002). These are, as best as we 
presently know, the physical parameters within which the 
asymmetry index is sensitive. 

These CDM merger fraction predictions, agree fairly 
well with our computed merger fractions at the highest 
redshifts. At the highest masses (M > 10 10 M© ) and 
brightest magnitudes (Mb < — 21) the difference between 
the observed merger fractions are within la of the Ben- 
son et al. (2002) model results at z > 2. Discrepancy is 
however present at some redshifts, and at various stellar 
mass and magnitude limits, particularly for the faintest 
and lowest mass systems. This discrepancy is highest at 
z ~ 1 for systems with low masses (M* < 10 95 M Q ) and 
faint magnitudes (Mb > —20), where we find differences at 
significances of > 3 — 4<r. At other redshifts we find that 
the differences are lower with significances of <~ — 2cr. 
These differences result from the fact that in hierarchical 
formation models, low-mass objects form by merging first, 
while it appears from our observations that lower mass 
galaxies tend to form continuously throughout the history 
of the universe, with a significant fraction of low-mass sys- 
tems not undergoing mergers at any redshift. 

We can also use merger fraction predictions from Khoch- 
far & Burkert (2001) to test qualitatively how well CDM 
theory predicts observed merger fractions as a function of 
redshift, stellar mass, and absolute magnitude. Khoch- 
far & Burkert (2001) suggest that at brighter limits the 
merger fraction power law fits have changing fo and toa 
values such that at higher mass and brighter magnitude 
limits fo goes down while toa goes up. In other words, the 
present-day merger fraction (fo) is lower for more massive 
galaxies, while the redshift evolution of the merger rate 
(771a) is steeper. We see this from the fits to the data 
(Tables 7 & 8) and also in terms of the mass evolution of 
the galaxies we study (§4.4). For the most massive and 
brightest systems we find that fo ~ and hia ^4 — 6, 
consistent with the prediction that these massive systems 
formed by merging early (Khochfar & Burkert 2001). 

The most interesting discrepancy between the computed 
merger fractions and those predicted by CDM is at a red- 
shift of z ~ 1. Although some of the significance in the 
difference (~ 4(7 ) can probably be accounted for by cos- 
mic variance, we speculate that it is unlikely to account 
for such a high z ~ 1 galaxy merger fraction. While it is 
possible that there exists a significant number of very low- 
surface brightness galaxies at z ~ 1 which we are missing, 
we are not likely missing normal galaxies at these redshifts, 
as these are clearly seen when lower redshift galaxies are 
placed at z <~ 1 — 2.5 (see §2.3.5). We therefore conclude 
that at z ~ 1 a physical effect is occurring which increases 
the merger fraction, or we are missing a population of 
galaxies for some unknown reason. The fact that CDM 
simulations, which have ad hoc prescriptions for creating 
star formation, and our merger fractions based on a small 
area of the sky, agree so well is suggestive that at least the 
most massive galaxies do indeed form by mergers. 

5. SUMMARY 

In this paper we present the first direct evidence for the 
hierarchical assembly of massive galaxies, as well as mea- 
surements of merger and merger mass fractions and rates 
at redshifts z > 1. We are able to measure this evolu- 



tion by identifying galaxies undergoing major mergers out 
to z ~ 3 using the asymmetry parameter as described in 
Conselice et al. (2000a) which is part of the CAS morpho- 
logical system (Conselice 2003). After simulating how our 
asymmetry measurements change due to decreased S/N 
and reduced resolution inherent at higher redshifts, inves- 
tigating completeness, the Eddington bias, and correcting 
for these effects, we are able to conclude the following: 

i) The merger history of field galaxies changes as a func- 
tion of absolute rest-frame magnitude Mb, stellar mass 
(M*) and redshift (z). For galaxies with Mb > —20 or 
M* < 10 95 M the merger fraction peaks at a value ~ 0.2 
near z ~ 1 and slightly declines at higher redshift. Fitting 
these lower mass and fainter merger fractions to a simple 
power-law of the form f = fo x (1 + z) mA we find power-law 
slopes, 777a ^2.5 — 3 out to z ~ 1 and toa ^0.5 — 1 from 
z ~ — 3. The corrected sizes of galaxies also become 
larger at lower redshifts (Figure 11) (see also Papovich et 
al. 2003). 

ii) We see a clear bifurcation in the merger fraction evolu- 
tion for the most massive, M* > 10 10 M , and brightest, 
Mb < —21, galaxies such that the merger fraction con- 
tinues to increase for these systems at higher redshifts, 
with peak values near 0.5 at z ~ 2.5 and low fractions 
~ at z ~ 0. Fits to the power-law evolution of the 
merger fractions for these galaxies reveals very steep slopes 
with m,A ~ 4 — 6. That is, luminous systems underwent 
more frequent major mergers at high redshifts than lower- 
luminosity galaxies, or their low-redshift bright and mas- 
sive counterparts. 

iii) By using results from N-body simulations of galaxies 
involved in major mergers we are able to convert merger 
fractions into co-moving volume merger rates, 4>m, finding 
a peak merger rate at z ~ 1 of <j>M ~ 2 x 10 6 mergers 
Gpc -3 Gyr -1 for systems with M B < -19. 

iv) Using stellar masses (M+) of the HDF galaxies mea- 
sured by Papovich (2002) we determined the evolution in 
the mass accretion rate per co-moving volume, the mass 
accretion rate per galaxy, and the fraction of mass and lu- 
minosity in galaxies undergoing major mergers as a func- 
tion of redshift. We find that the peak mass accretion 
rate per co-moving volume is ~ 6 x 10 6 M Gpc -3 year -1 
at z ~ 1. We also find that the mass accretion rate per 
galaxy, Mq increases as a function of redshift for all galax- 
ies, as does the fraction of mass in galaxies undergoing 
mergers. For galaxies with M* > 10 9 M Q we find a max- 
imum mass per galaxy accretion rate of Mq ~ 5.5 x 10 8 
M Q Gyr -1 at z ~ 2.5, with a rate evolution given by 
M G = 1.6 x 10 8 M Galaxy -1 Gyr" 1 (1 + z )0-99±o.32_ 
The fraction of galaxy stellar mass density involved in 
mergers also increases as a function of redshift, but much 
more rapidly, and with a higher maximum fraction, for the 
brightest and most massive systems. For galaxies with 
M B < -21 or M* > 10 10 M Q , the fraction of mass in- 
volved in mergers is ~ 0.5 at z ~ 2.5, demonstrating that 
at least half of mass in the most massive galaxies in the 
nearby universe were involved in major mergers ~ 9 Gyrs 
ago. 

v) Qualitative and quantitative comparisons of merger 
fractions with results from Cold Dark Matter (CDM) simu- 
lations are in relatively good agreement for the most mas- 
sive systems at z ~ 2.5. There is some discrepancy be- 



tween the models and the observed merger fractions for 
galaxies at low-masses and faint magnitudes, especially at 
z ~ 1. Our results are in agreement with the core re- 
sult of CDM structure-formation models in which massive 
galaxies form and evolve by merging. 

Changes in the merger history through time can also 
explain a host of galaxy phenomenon which we have not 
considered or discussed in this paper, including: variations 
in cosmic star formation history (e.g. Madau et al. 1998), 
the peak in density of active galactic nuclei (Boyle & Ter- 
levich 1998), the formation of black holes (Menou, Haiman 
& Narayanan 2001), and the evolution of supernovae and 
gamma-ray bursts. In general, our results are in agree- 
ment with the idea that major mergers have occurred in 
large numbers in the past, and that a large fraction of 
the most massive galaxies in the universe have formed by 
the merging of lower mass systems. While our results are 
for only a small area of the sky, future observations with 
the Advanced Camera for Surveys on the Hubble Space 
Telescope, such as the GOODS fields, will allow us to put 
firmer constraints on the merging history of galaxies, in- 
cluding determining how the merger rate varies as a func- 
tion of environment. 
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TABLE 1 
Simulation Results 



Redshift (z) 


Number Detected (% Detected) 


6(A)* 


<5(r) (arcsec) a 


0.5 


38 (1.00) 








1.0 


38 (1.00) 


0.05±0.09 


0.12±0.14 


1.5 


38 (1.00) 


0.04±0.07 


0.11±0.10 


2.0 


38 (1.00) 


0.04±0.08 


0.14±0.14 


2.5 


38 (1.00) 


0.06±0.09 


0.20±0.17 


3.0 


34 (0.89) 


0.06±0.12 


0.26±0.27 


4.0 


23 (0.61) 


0.08±0.13 


0.41±0.45 


5.0 


12 (0.32) 


0.11±0.16 


0.62±0.68 



a <5(A) and 5(r) are the differences between the originally measured asymmetries and radii for the 38 galaxies at 0.4 < z < 0.7 simulated at higher z and the 
resulting values measured when the galaxy is at the given higher Redshift (z). 
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TABLE 2 

Average Random Asymmetry Errors at Various Magnitude and Redshift Limits 



Redshift Range 


-19 < M B < -18 


-20 < M B < -19 


-21 < M B < -20 


M B < -21 


- 1.0 


0.02 


0.01 


0.01 


0.00 


1 - 2.0 


0.14 


0.07 


0.04 


0.03 


2 - 3.0 


0.18 


0.13 


0.08 


0.04 


3 - 3.5 


0.19 


0.17 


0.15 


0.10 
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TABLE 3 

Major Merger Candidates in the Hubble Deep Field North & 



NIC ID b 


R.A. (J2000) 


Dec. (J2000) 


Redshift 


A(B) rest d 


M B 


(B-V) rest 



1022 


12 


36 


:50 


2 


62 


■12 


:39. 


8 





.47 





40±0 


00 


-20 


.62 


0.50 


1335 


12 


36 


:41 


4 


62 


11 


:42 


5 





.55 





,75±0 


01 


-18 


.70 


0.71 


826 


12 


36 


39 


7 


62 


12 


:29 


.4 





.83 





,50±0 


.05 


-18 


.48 


0.19 


330 


12 


36 


:50 


.0 


62 


:13 


:51 








.85 





,64±0 


03 


-19 


93 


0.28 


40 


12 


36 


48 


.6 


62 


14 


:23 


2 





.95 





61±0 


02 


-19 


33 


0.31 


909 


12 


36 


:57 


.7 


62 


:13 


:15 


2 





.95 





,44±0 


02 


-20 


.46 


0.33 


488 


12 


36 


■48 


.6 


62 


:13 


28 


3 





.96 





,37±0 


20 


-20 


.36 


0.32 


1076 


12 


36 


:40 


.8 


62 


12 


:03 


1 


1 


.01 


o 


39±0 


.03 


-20 


23 


0.29 


152 


12 


36 


■48 


.3 


62 


:14 


:12 


4 


1 


.02 





,37±0 


.02 


-19 


60 


0.35 


1447 


12 


36 


:44 


.5 


62 


11 


:41 


6 


1 


.02 





,37±0 


.02 


-20 


.39 


0.33 


1552 


12 


37 


■01 


.6 


62 


:12 


:26 


8 


1 


.05 





,67±0 


.01 


-18 


.17 


0.59 


576 


12 


36 


49 


.1 


62 


:13 


:21 


9 


1 


.08 





,59±0 


.11 


-19 


.57 


0.64 


519 


12 


36 


■42 


.5 


62 


:13 


:0B 


.2 


1 


.10 





,51±0 


.09 


-19 


.37 


0.39 


1652 


12 


36 


49 


.7 


62 


11 


:49 





1 


.22 





71±0 


.01 


-18 


.25 


0.33 


1090 


12 


36 


56 


.6 


62 


:12 


:52 


.7 


1 


.23 





,54±0 


.02 


-19 


70 


0.36 


1568 


12 


36 


:54 


.3 


62 


12 


:02 


6 


1 


.23 





,36±0 


02 


-19 


.67 


0.54 


1141 


12 


36 


;50 


.0 


62 


■12 


:26 


3 


1 


.23 





,86±0 


04 


-19 


49 


0.37 


884 


12 


36 


.55 


.2 


62 


:13 


:09 





1 


.27 





39±0 


.02 


-19 


.25 


0.56 


577 


12 


36 


42 


.7 


62 


:13 


:06 





1 


.29 





41±0 


.16 


-20 


.58 


0.64 


872 


12 


36 


:51 


.6 


62 


:13 


:00 


3 


1 


.34 





68±0 


.03 


-18 


.56 


0.25 


360 


12 


36 


:52 


.7 


62 


:13 


:54 


8 


1 


.36 





,59±0 


05 


-22 


.17 


0.27 


368 


12 


36 


■47 


4 


62 


:13 


:37. 


.2 


1 


.41 





,45±0 


.04 


-18 


.09 


0.44 


1431 


12 


36 


:42 


.7 


62 


11 


:40 


8 


1 


.44 





,37±0 


12 


-18 


.31 


0.39 


1021 


12 


36 


:46 


.2 


62 


12 


:28 


.5 


1 


.45 





63±0 


.05 


-19 


.94 


0.36 


1520 


12 


37 


:00 


.8 


62 


12 


:27 


6 


1 


.50 





.43±0 


17 


-18 


.12 


0.39 


1188 


12 


36 


.54 


.0 


62 


12 


.35. 


2 


1 


.50 





,38±0 


.05 


-18 


20 


0.09 


1570 


12 


36 


:58 


.1 


62 


12 


:14 


3 


1 


.51 





39±0 


.07 


-18 


27 


0.33 


1512 


12 


36 


:52 


1 


62 


12 


:01 


2 


1 


.53 





,36±0 


09 


-19 


96 


0.66 


970 


12 


36 


:45 


.6 


62 


12 


:33. 


5 


1 


.56 





,65±0 


17 


-18 


49 


0.39 


1020 


12 


36 


46 


.2 


62 


12 


:29 


1 


1 


.75 





,84±0 


.06 


-19 


54 


0.21 


422 


12 


36 


43 


2 


62 


:13 


:19 


1 


1 


.78 





36±0 


07 


-18 


.33 


0.33 


957 


12 


36 


45 


.4 


62 


12 


:33. 


7 


1 


.81 





,49±0 


14 


-20 


12 


0.51 


1140 


12 


36 


49 


.9 


62 


12 


:27 


.0 


1 


.86 





36±0 


.22 


-18 


.92 


0.35 


732 


12 


36 


:44 


.4 


62 


12 


:44 


1 


1 


.98 





,40±0 


.11 


-21 


26 


0.35 


229 


12 


36 


:51 


3 


62 


14 


:11 


3 


2 


.01 





,46±0 


15 


-20 


80 


0.23 


466 


12 


36 


:53 





62 


:13 


:44 


2 


2 


.02 





66±0 


.15 


-21 


.21 


0.29 


1583 


12 


36 


54 


.8 


62 


12 


03 


2 


2 


.49 





,55±0 


.12 


-20 


.36 


-0.26 


1252 


12 


36 


:51 


.7 


62 


12 


:21 


4 


2 


.71 





63±0 


03 


-21 


.06 


0.38 


604 


12 


36 


■44 


.8 


62 


:13 


:07 





2 


78 


1 


,06±0 


10 


-21 


.35 




1358 


12 


36 


:45 


.3 


62 


11 


:52 


2 


2 


.80 





,45±0 


.05 


-21 


.97 


0.44 


522 


12 


36 


44 


.1 


62 


:13 


:10 


8 


2 


93 





36±0 


.04 


-22 


61 




813 


12 


36 


■47 


9 


62 


12 


:55. 


4 


2 


93 





,42±0 


.03 


-21 


.87 




1541 


12 


36 


48 


3 


62 


11 


:45 


8 


2 


.98 





,38±0 


.07 


-21 


17 





a These are the galaxies with rest-frame B-band asymmetries > 0.35 after correcting for systematic effects based on our simulations (§2.3.5). We also list only 
those systems at z > 2 that are bright than Mg — —20, as the identify of a major merger for fainter galaxies at z > 2 becomes less certain, and the implied merger 
fraction is more of a statistical quantity. 

b NIC ID is the ID number in the Dickinson et al. (2003b) catalog. 

c Rcdshifts are a mix of spectroscopic, when available, otherwise photometric. 



Rest-frame B-band asymmetries which have been corrected for systematic redshift effects (§2.3.5) 



18 



TABLE 4 

Average Rest-Frame B-band Asymmetries at Various Magnitude Limits 



Redshift Range 


M B < -IS 


M B < -19 


M B < -20 


M B < -21 


M B < -22 


- 1.0 


0.19±0.12 


0.21±0.12 


0.20±0.10 


0.16±0.06 


0.13±0.00 


1 - 2.0 


0.20±0.18 


0.21±0.16 


0.20±0.12 


0.23±0.13 


0.33±0.25 


2 - 3.0 


0.15±0.19 


0.17±0.17 


0.23±0.18 


0.35±0.22 


0.28±0.07 


3 - 3.5 a 


0.14±0.19 


0.13±0.19 


0.16±0.18 


0.17±0.17 


0.33±0.19 



a At the restshift range 3 - 3.5 we are no longer sampling the rest-frame B band morphologies of these galaxies but are viewing them in their near UV. 
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TABLE 5 

Inferred Merger Fractions (F merg er) as a Function of Luminosity and Redshift 8 



Redshift 
M B = -18 


Fmerger 
-18 


Num. Gal. 
-18 


Redshift 
-19 


Fmerger 
-19 


Num. Gal. 
-19 


Redshift 
-20 


Fmerger 

-20 


Num. Gal. 
-20 


Redshift 
-21 


Fmerger 

-21 


Num. Gal. 
-21 


0.58 


0.04 


51 


0.57 


0.04 


26 


0.59 


0.07 


15 


0.65 


0.00 


4 


1.10 


0.14 


142 


1.10 


0.18 


87 


1.10 


0.17 


41 


1.10 


0.08 


12 


1.73 


0.14 


93 


1.75 


0.09 


53 


1.84 


0.10 


20 


1.83 


0.17 


6 


2.41 


0.09 


183 


2.46 


0.10 


116 


2.46 


0.18 


56 


2.55 


0.4 


20 



a For each limiting magnitude, printed below each quantity, the fraction of galaxies consistent with undergoing a major merger is listed (F mer g er ), as is the total 
number of galaxies within the redshift and magnitude limit range (Num. Gal.). The redshifts listed arc the averages within the limits: 0.4 < z < 0.7,0.7 < z < 
1.3, 1.3 < z < 2, 2 < z < 3. 
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TABLE 6 

Inferred Merger Fractions (F merg er) as a Function of Limiting Mass and Redshift 3 - 



Redshift 


Fmerger 
8 


Num. Gal. 


Redshift 


Fmerger 


Num. Gal. 


Redshift 


Fmerger 


Num. Gal. 


Redshift 


Fmerger 


Num. Gal. 


log (M.) = 8 


8 


9 


9 


9 


9.5 


9.5 


9.5 


10 


10 


10 


0.59 


0.06 


117 


0.58 


0.04 


46 


0.57 


0.07 


27 


0.60 


0.06 


16 


1.10 


0.13 


260 


1.10 


0.14 


111 


1.10 


0.15 


61 


1.00 


0.06 


32 


1.71 


0.14 


118 


1.73 


0.17 


77 


1.77 


0.15 


26 


1.81 


0.11 


9 


2.37 


0.09 


150 


2.40 


0.13 


101 


2.43 


0.18 


51 


2.56 


0.5 


10 



a For each limiting lower mass the number fraction of galaxies consistent with undergoing a merger (F merger ) is listed as is the total number of galaxies within the 
redshift and magnitude limit range (Num. Gal.). The redshifts listed are the averages within the limits: 0.4 < z < 0.8, 0.8 < z < 1.4, 1.4 <z<2,2<z<3. 
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TABLE 7 

Merger Fraction Fits Using Absolute Magnitude Limits 3, 



Magnitude Limit (Mb) 


Using A 


Using A 


2^1 


z ~ 2 


constant fo 


constant fo 




fo 


m A 


m A 


m A 


fo 


m A 


-18 


0.06±0.05 


0.5±0.6 


2.5±0.3 


2.0±0.1 


0.021 


1.5±0.2 


-19 


0.07±0.08 


0.4±0.8 


2.8±0.3 


1.9±0.3 


0.021 


1.5±0.2 


-20 


0.07±0.05 


0.7±0.6 


2.8±0.1 


1.8±0.3 


0.021 


1.8±0.1 


-21 


0.004±0.001 


3.7±0.3 


1.4±1.B 


1.9±0.1 


0.021 


2.3±0.1 



a This table shows the merger fraction fitted parameters using the listed lower absolute magnitude limits (M B ). The merger fractions are fitted in four different 
ways using f m (A,Mg, z) — fo X (1 + z) rn A . The fitted form using only the fractions computed using the asymmetries out to z ~ 3 is listed in the 'Using A' column. 
The fitted merger fraction slope out to z ~ 1 and 2^2 are also listed. The final two columns shows merger fraction fits when using the 2^0 merger fraction 

point from Patton et al. (1997) and holding fg to match this zero redshift point. 
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TABLE 8 

Merger Fraction Fits Using Stellar Mass Limits' 1 



Mass Limit (log(M*)) 


Using A 


Using A 


z ~ 1 


z ~ 2 


constant fo 


constant fo 




to 


m A 


m A 


m A 


fo 


m A 


8.0 


0.07±0.03 


0.4±0.6 


2.4±0.1 


2.0±0.1 


0.021 


1.5±0.2 


9.0 


0.06±0.03 


0.9±0.7 


2.5±0.3 


2.2±0.1 


0.021 


1.7±0.1 


9.5 


0.06±0.01 


0.9±0.3 


2.6±0.0 


2.1±0.1 


0.021 


1.9±0.1 


10.0 


~ 


5.9±1.3 


1.7±0.3 


1.6±0.1 


0.021 


2.4±0.1 



a This table shows the merger fraction fitted parameters using the listed lower stellar mass limits (log(M*). The merger fractions are fit in four different ways 
using f m (A,M*, z) — f'o X (1 + z) rn A . The fitted form using only the fractions computed using the asymmetries out to z ~ 3 is listed at the 'Using A' columns. The 
fitted merger fraction slope out to z ~ 1 and z ~ 2 are also listed. The final two columns shows merger fraction fits when using the 2^0 merger fraction point 

from Patton et al. (1997) and holding fg to match this zero redshift point. 
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TABLE 9 

Merger and Mass Accretion Rates and Densities' 1 



Redshift Range (f> M (Mergers Gpc~ 3 Gyr" 1 ) 13 p M (M Gpc- 3 ) c M (Mq Gpc~ 3 yr~ 1 ) d M G (M Gyr- 1 ) 

0.4-0.8 3.9x10 s 1.6X10 15 1.8x10 s 2.0x10 s 

0.8-1.4 2.2x10 s 5.4X10 15 6.0x10 s 4.2x10 s 

1.4-2.0 5.4x10 s 1.8X10 15 2.0x10 s 3.7x10 s 

2.0-3.0 7.3x10 s 3.1X10 15 3.5x10 s 5.5x10 s 



a These merger and mass accretion rates and densities are computed for systems with Mg < — 19 and assuming that a major merger can be identified though the 
asymmetry technique for 900 Myrs. For the mass accretion rate we assume that all galaxies identified as major mergers had initial masses with a 1:1 ratio. 
b Merger Rate 0m 
c Mass Accretion Density pjj 
d Mass Accretion Rate Density M 
°Mass Accretion Rate Galaxy -1 M(-; 
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1 2 3 

Redshift (z) 

Fig. 1. — Plot of absolute magnitude vs. redshift for galaxies in the Hubble Deep Field North. The red symbols are for objects with 
(B— V)rcst > 0.5 and blue symbols are for (B— V) rcs t < 0.5. The large symbols are objects with confirmed spectroscopic redshifts while 
the smaller symbols are those objects with photometric redshifts. The vertical line at z ~ 2.5 is the limit we use for obtaining reliable 
morphological parameters and galaxy detections, and the solid and dashed horizontal lines shows the Mb < —18 and Mb < —19 lower 
luminosity limits we use throughout this paper (see text). 
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Fig. 2. — Plot of the asymmetries computed in the I band for the 38 galaxies in the HDF with Mb < —18 and at 0.4 < z < 0.7 at various 
fractions of the Kron multiplier (k r ) times semi-major axis (a) radii. The solid box is the average asymmetry value at every radius, while 
errorbars represent the 1 a variation of these averages. 
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Fig. 3. — The distribution of asymmetries for all 1212 galaxies observed in the HDF as measured in the F814W (I) band and the F160W 
(H) band after using five different background positions for computing the background correction. 
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Fig. 4. — Results of simulating the 38 HDF galaxies at 0.4 < z < 0.7 and Mb < —18 to higher rcdshifts and then redoing the entire 
asymmetry analysis from detection to measurement. The left panel shows the average and la variation of the asymmetry difference between 
the original asymmetry and the simulated asymmetry (A org - A s i m ). Also plotted on the right side of the left panel is the fraction of the 38 
galaxies that are detected through SExtractor after performing the simulation. The right panel shows the average and la variations of the 
difference, in arcseconds, between the original radii measured (Rad org ) and the simulated radii (Rad s i m ). 
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Fig. 5. — Plot of the rest frame B-band asymmetry measurements A(B re st) as a function of redshift (z). The size of the plotted point is 
proportional to (left) the luminosity (magnitude) of the galaxy and (right) the stellar mass of the galaxy. The smallest points are objects 
with a. -18 > M B > -19, b. 8 < M» < 9. The larger symbols are for galaxies with a. -19 > M B > -20, b. 9 < M* < 9.5, and a. -20 > 
M B > -21, b. 9.5 < M» < 10. The largest symbols are for galaxies with a. -21 > M B > -23 and b. M* > 10. The solid line is A(B Iest ) 
= and the dashed line is the limit for mergers at A(B m ergcr) = 0.35. The solid diagonal lines originating at A(B re st) = and A(B mer ger) 
= 0.35 show the average random error on the asymmetry index as a function of redshift and limiting magnitude, such that each line is either 
A(z) = 0.35 + error(z) or A(z) = 0— error(z). These lines are for galaxies at magnitudes, from the nearest to their respecitve horizontal line 
outward, M B < -21, -21 < M B < -20, -20 < M B < -19, and -19 < M B < -18. 
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Fig. 6. — Distribution of asymmetries as a function of absolute magnitude and rcdshift. The top four numbers are, clockwise top-left, the 
mean asymmetry, its 1 <r variation, its skewness and kurtosis values. The bottom left number is the difference in the average random error 
for the galaxies in that particular magnitude and redshift bin and the previous lower redshift bin at the same magnitude. The bottom right 
number is the significance, based on Monte Carlo simulations discussed in §3.2, that an increase in asymmetries are not due to an increase in 
random errors. The dashed vertical line shows the A mcrKor limit. 
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Fig. 7. — Images of all galaxies in the Hubble Deep Field North with Mb < —20, between < z < 1, ordered by increasing asymmetry, as 
seen in F814W. The upper number in each panel is the Mb of each galaxy and the bottom number is its redshift, while the number on the 
right hand side is the galaxy's A(B) value. Galaxies consistent with being mergers based on their asymmetries have a solid dark box in the 
right hand corner of their respective panels. 



continued - same as part a. expect all galaxies with Mb < —20, between 1 < z < 2 arc shown. 
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Fig. 7. — continued - 



same as part a. expect all galaxies with Mb < —20, between 2 < z < 3 are shown. 
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Fig. 8. — Plot of the average rest frame B-band asymmetry for galaxies at different magnitudes limits as a function of rcdshift. Variations 
of these averages are listed in Table 4. 
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Fig. 9. — The (B— V) color distribution of galaxies plotted as a function of redshift (z). The larger boxes are galaxies at Mb < —20 while 
the small dots are galaxies with Mb > —20. The solid lines show synthesis stellar population colors with ages of 0.1, 1 and 3 Gyr for systems 
that formed in bursts. 
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Fig. 10. — Relationship between the rest frame B-band asymmetry A(B rcB t) and absolute magnitude (Mb) for HDF galaxies at two different 
redshift bins: z < 1.5 and 1.5 < z < 3. The color of each point is dictated by the rest frame (B— V) color of each galaxy, such that blue points 
are systems with (B— V) < 0.5 and red points are for galaxies with (B— V) > 0.5. 
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Fig. 11. — The distribution of HDF galaxy radii (0.5 X k r x a), in the rest-frame B-band to z < 2.5 and observed H-band at z > 2.5, 
corrected for redshift effects (§2.3.5 & Table 1) plotted as a function of redshift (z). The size of the symbol is proportional to luminosity, such 
that the smallest symbols are for galaxies with —20 < Mb < —18 and the two larger symbols are for systems with —22 < Mb < —20 and 
M B < -22. 
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Fig. 12. — The rest-frame (B— V) color verses asymmetry diagram plotted at four different luminosity intervals listed at the top of each 
panel. Each panel also plots galaxies divided into different redshift ranges at each luminosity. The diagonal dashed line is the relationship 
between (B-V) color and the rest-frame B-band asymmetry yr(B rC st) found for nearby normal galaxies (Conselice et al. 2000a). This line is 
extrapolated to bluer colors by the solid diagonal line. The vertical long-dashed line shows the ^4 m crgcr — ^(B-V = 0.4) + 3<t(B-V) limit 
we use for identifying mergers (see text). 
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Fig. 13. — Plots of the evolution of merger fractions as a function of redshift. The solid large round symbols are merger fractions computed 
using the asymmetry technique. The green crosses are merger fractions computed after using the +lcr and -1<7 values for each galaxy's 
asymmetry. The other symbols are defined in the legend printed on each panel, and are merger fractions found from galaxy pairs by Patton et 
al. (1997), LeFevre et al. (2000) and Carlberg et al. (2000). The straight blue line shows the merger fraction fit in the form f = f o X (l + z) mA 
when only using merger fractions computed by using the asymmetry derived fractions, and the cyan line is this fit when using the asymmetry 
merger fractions and holding fo = 0.021, the z ~ value found by Patton et al. (1997). The solid red line and points are merger fractions 
selected in an analogous way from the semi-analytic CDM simulation results of Benson et al. (2002). 
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Fig. 14. — An analogous plot of Figure 13 except that the merger fractions are selected based on stellar mass limits, rather than absolute 
magnitudes. 
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Fig. 16. — The stellar mass accretion rate per galaxy Mq, in units of Mq Galaxy 1 Gyr 1 plotted as a function of rcdshift (2) for gal 
with stellar masses greater than 10 9 , 10 9 ' 5 and 10 10 M . The dashed lines arc fits to the accretion rates listed in equations 1-3. 
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Fig. 17. — Plots of the fraction of stellar mass and galaxy luminosity involved in mergers as a function of redshift and magnitude (left 
panels) and stellar mass (right panels). The different lines show the fraction of mass and rest-frame B-band luminosity involved in mergers 
at different luminosity and stellar mass upper limits. 



